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We study the mechanics of nematically ordered semiflexible networks showing that they, like 
isotropic networks, undergo an affine to non-affine cross-over controlled by the ratio of the filament 
length to the non-affinity length. Deep in the non-affine regime, however, these anisotropic networks 
exhibit a much more complex mechanical response characterized by a vanishing linear response 
regime for highly ordered networks and a generically more complex dependence of the shear modulus 
upon the direction of shear relative to the nematic director. We show that these features can be 
understood in terms of a generalized floppy modes analysis of the non-affine mechanics and a type 
of cooperative Euler buckling. 

PACS numbers: 62.20.D- 87.16.Ln 82.35.Lr 



The mechanical properties of living eukaryotic cells are 
controlled by a low-density cross-linked network of semi- 
flexible protein filaments pQ. These filaments, composed 
primarily of F-actin, are densely cross-linked on the scale 
of their thermal persistence length. Understanding the 
material properties of such semiflexible gels is not only 
a forefront problem in biophysics but also presents a 
broader challenge in polymer physics. These cytoskeletal 
semiflexible networks can support elastic stress in both 
filament stretching and bending. This is in contradis- 
tinction to typical polymeric gels for which rubber elas- 
ticity theory applies [2]. Those flexible polymer systems, 
where the thermal persistence length of the constituent 
filaments is significantly smaller than the mean distance 
between cross-links along a given filament, store elastic 
energy only in the stretching of filament random walks 
between consecutive cross- links. 

Previous work [3j [4] has shown that, as a function 
of increasing network density or filament bending stiff- 
ness, semiflexible networks admit a sharp cross-over from 
a bending dominated elastically compliant regime that 
sets in for networks at densities above the rigidity per- 
colation transition [5 , to stiffer networks where elastic 
strain energy is stored primarily in the stretching of fila- 
ments. Moreover, under uniformly applied strain at the 
boundaries, the geometry of the deformation field in the 
softer, bending dominated regime is spatially heteroge- 
neous (i.e., non-affine) and shows large deviations from 
the affine deformation prediction of continuum elasticity 
theory [6 j. This cross-over is controlled by the ratio of 
the filament length L to the so-called non-affinity length 
A. 

This earlier work concentrated on the mechanics of 
statistically isotropic filament networks. Cytoskeletal 
networks however, are dominated by oriented stress 
fibers and typically have lower symmetry. Motivated 
by this point we explore in this Letter the mechanics 
of anisotropic networks by considering random networks 



with a nonuniform probability distribution of filament 
orientations. Although these anisotropic networks, char- 
acterized by both a filament density and nematic or- 
der parameter, undergo a non-affine to affine cross-over 
that is also controlled by L/\, we find that the mechan- 
ics of anisotropic semiflexible networks are considerably 
more complex. Specifically, we find that: (i) their lin- 
ear response regime narrows dramatically at high ne- 
matic order, and (ii) their shear response in the non- 
affine regime develops a complex dependence upon the 
angle between the shearing direction and the nematic di- 
rection that cannot be understood in terms of continuum 
linear elasticity. We can account for both of these re- 
sults using a self-consistent model of cooperative Euler 
buckling closely related to the "floppy modes" analysis 
of non-affine mechanics by Heussinger and Frey [4 . The 
analysis of anisotropic semiflexible networks more gen- 
erally provides deeper insight into the mechanics of the 
non-affine regime. In spite of the greater apparent com- 
plexity of anisotropic networks, their order makes them 
more amenable to analysis than the isotropic networks 
where the affine /non-affine cross-over was discovered. 

We began by numerically studying two-dimensional fil- 
ament networks composed of identical filaments of length 
L and having bending and stretching moduli n and fi re- 
spectively. To produce a system of N filaments with sta- 
tistically uniform density and filament anisotropy char- 
acterized by nematic order parameter 5, one endpoint of 
each filament was laid down at random in a square box of 
edge length W and area W 2 = A, and assigned an orien- 
tation chosen from the distribution P(0) with respect to 
nematic direction h. Crossing points between filaments 
were treated as incompliant crosslinks. Filaments were 
added until the system reached the selected density. A 
representative example of such a network, a magnified re- 
gion thereof, and the distribution of filament angles are 
shown in Fig. 0. The filaments were discretized such 
that nodes were placed at all crosslinks and at regular 
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FIG. 1: Top left: A typical nematic network with diagram 
defining the shear orientation angle (j). Bottom left: The an- 
gular distribution P(9) of filament orientations relative to the 
nematic director. Right: A closeup of a network (shear direc- 
tion indicated by arrows) showing the partitioning of elastic 
energy into filament stretching (red) and bending (blue). 



intervals between crosslinks to allow for bending. The 
energy of the network is given by 



n c (0)=2pL 2 / dip 8mi/>P(0 + ip). 



(2) 



Integrating over the orientation of that original filament, 
we find that the mean distance between cross-links is 
given by 



tc = L 
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where the angle brackets denote averaging over 6 with 
weight P{6). For an isotropic gel, Head et. al found that 
non-affine to affine cross-over is controlled by L/X where 
A = £ c (£c/£b) z with z = 1/3 [3 j. Scaling arguments [3] 
and mean field theories [8] suggest z = 2/5, 1/4 respec- 
tively. 

Two-dimensional elastic continua with nematic order 
may be characterized by four independent elastic con- 
stants [9] . Writing the elastic constant tensor in the usual 
way, E = \ j cPxCijkiUijUki [10 , we can identify these 
as Cini, C2222, C1212, and 61122- Assuming affine defor- 
mation, these quantities are all linearly proportional to 
the stretching modulus of the constituent filaments via 
the relation 
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where s is the length of a segment of rest length sq> /3 
is the angle between adjacent segments on one filament, 
and £0 is the average rest length of the two segments 
surrounding an angle spring. We used the Lees-Edwards 
method [7] to shear the system with periodic boundary 
conditions. 

The networks are characterized by a nematic order 
parameter S and a density, measured in terms of the 
mean distance between cross-links along a filament, £ c . 
The mechanics of the monodisperse filaments is set by 
£b = \JkI [i. Finally, all energy scales are set by a single 
filament elastic modulus; we measure all moduli in terms 
of p. 

For a given P(0), the nematic order parameter S = 
dO P(9) cos 26 describes the degree of anisotropy; S = 
corresponds to an isotropic system, and S = 1 is a sys- 
tem with all filaments aligned. We compute the mean 
distance between cross-links from the filament density 
p = N/A as follows: Consider a straight filament ly- 
ing at an angle 6 with respect to the nematic director. 
Given a second filament at some angle ip with respect to 
the first, the probability of their crossing is L 2 \ smip\/A. 
Integrating over ip with a weight P(0 + ip) we find that 
the mean number of filaments crossing the original one 
lying at angle with respect to n is Poisson distributed 
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where e(0) is the direction along a filament oriented at 
angle and n , the number of cross-links per filament, is 
summed over using their Poisson distribution p n (0). We 
focus exclusively on the shear modulus G(0), which, for 
a given shear orientation angle </>, the angle between the 
displacement and nematic directions as shown in Fig. [l] 
takes the form, using Eq. Q, of 

GaffmeW = G affin e(0) [l - 8 sin 2 (j) cos 2 +T sin 2 (j) cos 2 (j) 

(5) 

with 



T = upL / d<9P(<9) V 
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We refer to the above as the affine prediction for the shear 
modulus of our anisotropic solids. The angular depen- 
dence of G a ffine(</>) has a simple interpretation. A shear 
deformation with displacement along x is, to linear order, 
equivalent to stretching and compressing along the ±7r / 4 
directions with respect to x. Filaments oriented at an- 
gles ±7r/4 with respect to x will be stretched/compressed 
along their axes, and thus have the most strain energy. 
Thus networks sheared so that the displacement and ne- 
matic director are separated by these angles will have 
maximal shear moduli. 
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The depression of the shear modulus G below its amne 
value G affi ne serves as a mechanical measure of the de- 
parture from affinity. Geometric measures of this de- 
parture [3j [11] measure the spatial heterogeneity of the 
strain field as a function of length scale. Following Head 
et al. [3 we adopt as our geometric measure of nonafflnity 

A(r) = ((w — ^affine) 2 )(j) 7 wnere ^ is the strain-induced 
change in angle between an arbitrary axis and the line 
joining two network points separated by a distance r, 
^affine is the change of this angle given purely amne de- 
formations, and the average is taken over all point pairs 
with separation r in a network. 

Iff 1 I 1 I 1 I 1 I 1 I 1 F= 
I #-# S = 0.0, (J) = I 




FIG. 2: Normalized shear modulus G/G a m n e versus L/X for 
varying degrees of anisotropy S and shear orientation angle 
<f>. For highly anisotropic gels non-affine networks is there 
a large shear orientation angle dependence to the deviation 
from the amne shear modulus prediction. Inset: The geomet- 
ric nonamnity measure A as a function of length scale for six 
different cases, including the S = 0.7 network at various val- 
ues of 0, showing that A, unlike the modulus, is insensitive 
to the shearing angle. Moduli measured at 7 = 10 -3 . 

Fig. [2] plots the normalized shear modulus G/G a ffi ne as 
a function of L/X for various values of S and <ft. These 
data show reasonable collapse onto a universal curve, 
demonstrating that the amne to non-affine cross-over per- 
sists for anisotropic networks. The shear modulus of pre- 
sumably amne networks (L/X > 15) nearly saturates the 
amne prediction and has a shear orientation angle depen- 
dence as predicted by the amne theory. This can be seen 
seen by the collapse of the ratios G/G a ffi ne for various 
values of (j) in the amne regime. Deep in the non-affine 
regime, however, the measured shear modulus of highly 
anisotropic networks (S = 0.7) varies more strongly with 
<j) than can be understood by the amne theory. More- 
over, this additional angle dependence breaks the univer- 
sal data collapse observed in isotropic and more weakly 
anisotropic networks. This breakdown of continuum lin- 
ear elasticity is also seen in the divergence of A(r) as 
r — » (inset) [3]. This geometric measure is independent 
of the shearing angle. Thus, the degree of heterogeneity 



of the deformation field depends solely on L/X, but the 
mechanics of the anisotropic networks in the non-affine 
regime appears to have a novel and more complex depen- 
dence on <j). 
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FIG. 3: Shear modulus as a function of the nematic order 
parameter S (for fixed filament density) normalized by its 
value at S — and L/X ^ 5 . 

To understand the mechanics of anisotropic networks 
deep in the non-affine regime, we turn to an extension of 
the "floppy modes" model of Heussinger et. al [4 . This 
model assumes that the deformation energy is stored pre- 
dominantly in filament bending. For a filament segment 
of length i this energy is i^bend — ^na/^ 3 , where S na is 
the typical size of axial displacements. For larger values 
of i, it is energetically favorable for segments to bend, 
thus inducing a small amount of bending on connecting 
filaments; however, below some length £ m { n , the energetic 
cost of bending becomes too high. The displacement of 
those filaments results in stresses in their neighbors caus- 
ing them to bend. One self-consistently determines £ m i n 
by balancing the bending energy of these small segments 
with the total bending energy of the sea of longer fila- 
ments to which they are connected. This balance deter- 
mines the elastic energy storage and thus the modulus. 

We consider an idealized anisotropic gel with two 
classes of filaments: (i) nematic filaments oriented at 
= and (ii) impurity filaments oriented at = tt/2. 
The fraction of nematic filaments 25 — 1 is set to repro- 
duce the correct nematic order parameter. Now the bal- 
ance of energy of nematic filaments with their impurity 
neighbors and impurity energy with nematic neighbors 
results in two equations that must be solved simultane- 
ously: 

K() 2 C°° K() 2 

js^ = n c,i / d£jP(£j)^ and 

^min,N «^min,i ^1 

KC) 2 r°° KC) 2 

W^ = n c ,N de N p(e N )-p, (?) 

^min,I «^min,N N 
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where P(£) is the distribution of segment lengths nematic 
(N) and impurity (I) filaments. The solution, as in the 
isotropic case, gives G ~ p 7 , but also predicts the de- 
pendence of the shear modulus on S deep in the non- 
affine phase. In Fig. [3] the anisotropic floppy mode the- 
ory (sold line) agrees well with the simulation data (red 
dots), while the affine prediction (green dashed line) does 
not. 
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FIG. 4: (a) Difference in G for = ±7r/4 at 5 = .7 as a 
function of L/A. The shear modulus was measured strain 
7 = .001; (b) G/G a ffine vs. 7 for S = .7 for various values of 
L/A and 0; (c) Fraction of energy in stretching as a function 
of shear for the same networks. 

This result does not explain the unexpectedly stronger 
<fi dependence in the non-affine regime. To account for 
this effect, we must turn to the nonlinear response of 
anisotropic networks in the non-affine regime. Fig. |4ja) 
plots the difference between G measured at <j) = 7r/4 and 
(j) = — 7r/4 at 7 = 0.001 as a function of L/A for a network 
with S = 0.7. G 1 is computed as the second derivative of 
the energy density with respect to 7. As these two orien- 
tations are equivalent under a mirror reflection, according 
to continuum linear elasticity AG 7 should vanish. In- 
deed it does so as we reach the affine regime: L/A ^> 1. 
For non-affine networks, however, AG 7 becomes large. 
Examining Fig. [4^b) we see the origin of this particular 
deviation from continuum linear elasticity for non-affine 
gels (L/A = 7.9). AG 7 increases with 7, demonstrating 
that it is an inherently nonlinear effect. Specifically, the 
modulus for (j) = — 7r/4 (dotted black line) monotonically 
softens relative to that of <fi = 7r/4 (solid black line) as 
7 increases. For an affine network (L/A = 29) these two 
shear moduli (green solid and dotted lines) are equal for 
all strains explored showing that affine networks have a 
significantly larger linear response regime. 

The nonlinearity that leads to the difference AG 7 can 



be traced to cooperative Euler buckling of the network 
that leads to the nonlinear softening of the shear modu- 
lus at (j) — —tt/4 since, at this network orientation rela- 
tive to the shearing direction, more filaments are placed 
under compressive stress. Moreover, these filaments have 
larger values of £ c allowing for buckling at smaller strains. 
To demonstrate this, we plot in Fig. |4^c) the fraction of 
elastic energy stored in stretching as a function of both 
strain and L/A. In the non-affine regime L/A = 7.9, 
there is generically a smaller fraction of elastic energy 
stored in stretching, as expected from previous work, but 
the fraction of stretching energy decreases dramatically 
with strain as the network is sheared at angle <fi = —it /A 
(dotted black line) applying compressive stresses to the 
filaments along ft. Shearing along <\> = 7r/4 puts the ne- 
matically aligned filaments under tension. The filaments 
now under compression are generically cross-linked on a 
much finer scale so that the network now collectively re- 
sists buckling. 

To estimate the shear at which we expect AG 7 > 0, we 
note that Euler buckling of a filament of length £ occurs 
above a critical compressive load F c = 7r 2 tv/£ 2 [6 j. For 
an affine deformation, F = with the compressive 

strain along <p = — tt/4. We find a critical strain j c = 
2tt 2 (£i ) /£) 2 to cause a simply-supported segment of length 
£ to buckle. To find upper and lower bounds for 7 C , we 
imagine that buckling sets in on scales £ c < £ < L. For 
L/A = 7.9 (Fig. g we find that 3 x 10~ 5 < j c < 10 
reasonable agreement with the data. 
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We have found that anisotropic networks, much like 
isotropic networks, undergo an affine to non-affine cross- 
over that is controlled by L/A. The geometric measure 
of no n- affinity depends on L/A in a manner independent 
of anisotropy. The mechanics of the non-affine regime of 
anisotropic networks, however, is quite complex, charac- 
terized by a vanishing linear response regime [12] in the 
limit of increasing nematic order parameter. The over- 
all dependence of the network's shear modulus upon S 
can be understood in terms of a generalization of the 
floppy modes picture. The vanishing of the linear re- 
sponse regime reflects the appearance of a type of co- 
operative Euler buckling. We note that, while affine net- 
works of arbitrary anisotropy can be understood in terms 
of elastic continua, non-affine and anisotropic networks 
allow for a much more complex set of mechanics. This 
result suggests that the mechanics of ordered cytoskeletal 
structures must be modeled with particular care in the 
non-affine regime. Finally, we propose that insights into 
the highly nonlinear mechanics of anisotropic non-affine 
networks should lead toward a more complete theory of 
the mechanics of the non-affine regime, particularly in 
the nonlinear regime. 
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